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AC Stark effect in ThO H 3 A\ for the electron EDM search 
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A method and code for calculations of diatomic molecules in the external variable electromagnetic 
field have been developed. Code applied for calculation of systematics in the electron’s electric dipole 
moment search experiment on ThO H 3 Ai state related to geometric phases, including dependence 
on O-doublet, rotational level, and external static electric field. It is found that systematics decrease 
cubically with respect to the frequency of the rotating transverse component of the electric field. 
Calculation confirms that experiment on ThO H 3 Ai state is very robust against systematic errors 
related to geometric phases. 


The experimental measurement of a non-zero electron 
electric dipole moment (eEDM, d e ) would be a clear sig¬ 
nature of physics beyond the Standard model [INi. The 
current limit for eEDM, |d e | < 9 x 10 -29 e-cm (90% 
confidence), was set with a buffer-gas cooled molecular 
beamj^-Q of thorium monoxide (ThO) molecules in the 
metastable electronic H 3 Ai state. It was shown that 
due to existence of closely-spaced levels of opposite par¬ 
ity of D-doublet the experiment on ThO is very robust 
against a number of systematic effects related to mag¬ 
netic fields[§1, [9] or geometric phases[lo|]. However, the 
upper and lower D-doublet states have slightly differ¬ 
ent properties and systematic effects related to magnetic 
field imperfections and geometric phases can still mani¬ 
fest themselves as a false eEDM. The dependence of g- 
factors of the ThO H 3 Ai state on D-doublets and exter¬ 
nal electric field was considered in Ref. )9j. The aim of 
the present work is to consider geometric phase shifts. 

Following the computational scheme of 1, U|, the en¬ 
ergies of the rotational levels in the H 3 A\ electronic state 
of the 232 Th 16 0 molecule in external static electric £ = 
£z and magnetic B = Bz fields are obtained by numeri¬ 
cal diagonalization of the molecular Hamiltonian (H mo i) 
over the basis set of the electronic-rotational wavefunc- 
tions 'hfi0^ -in (a,/3). Here To is the electronic wave- 
function, 6 J MSl (a,P) = i/( 2J + l)/ 47 rD^ n (a,/ 3,7 = 0) 
is the rotational wavefunction, a, /3 ,7 are Euler angles, 
and M (f V) is the projection of the molecule angular 
momentum J on the lab z (internuclear h) axis. De¬ 
tailed feature of the Hamiltonian is described in ti|. In 
the paper the M = ±1 states which represent interest 
for eEDM search experiment are considered. For elec¬ 
tric field £ = 20 — 200 V/cm, used in the experiment, 
lower rotational levels with M ^ 0 can be labeled by 
| J,M,Q > quantum numbers. States |J, M=1,0=1 >, 

| J,M= — l,f]= — 1 > correspond to the upper and 
| J,M= — l,fi=l >, \J,M=l,Cl= — 1 > to the lower 
H-doublet levels. External magnetic field removes the 
degeneracy between D-doublet components: A E u = 
E(\J , M= 1, ft=l >) - E(\J, M= - 1, Q= - 1 >), A E l = 
E(\J, M=l, 0= — 1 >)-E(\J,M=-l,n=l >). The rel¬ 


evant energy levels can be seen in Figure 2 of Ref. B] or 
Figure 3 of Ref. 0. Provided g-factors for upper and 
lower D-doublet levels are close enough A E u and A E l 
remain equal unless both parity and time reversal sym¬ 
metries are violated. The difference in splitting gives the 
value for eEDM d e = ^ AE 4 £ AE ^ , here £ e s = 81.5GV/cm 
Bill is the effective internal electric field. However, 
there are systematic effects which can give additional en¬ 
ergy shifts 6AE 1 and 5AE U for A E l and A E u which 
manifest as a false eEDM. This leads to a systematic er¬ 
ror Sd e (sys) = 5AE 4 £ SAE ■ It is also useful to consider 

systematic effects 5d e (sys) = SAE * > related to one of 
the fl-doublet component . One of the effect is the in¬ 
teraction with transverse component of the electric field 
£{t) = £ j_(icos(o;_Lf) + ysin(oj±t)) which appears due to 
a spatial inhomogeneties in the applied electric field 0 ] . 
Let us consider this effect. 

The corresponding part of the Hamiltonian is 


H tf = -d ■ £(t) = -£ ± /2(d + e~ iu ^ + d-e iu ^), (1) 


where d± = d x ± id y . It is more convenient to describe 
the interaction of the molecule with the quantized elec¬ 
tromagnetic fields. The corresponding Hamiltonian is 


Hint = hui±a + a — 


27T?kUj 


V 


(d+a + + d-a), ( 2 ) 


where a + and a are creation and annihilation operators, 
V is a volume of the system. To work with Hamiltonian 
(ED one need to add the quantum number |n >, where 
n = jp —7 is number of photons. The approach is similar 
to the formalism outlined in Ref. 0 - For this paper 
we consider the case £± « £, such that the additional 
energy shifts can be calculated in the framework of the 
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second order perturbation theory. 

6AE U = £i/4x 

| < J,M=l,n=l\d + \J',M=Q,n' > | 2 
E(J, M=l, n=l) - E(J', M= 0, Q/) + huj ± 


E 


| < J,M=l,fl=l\d_\J',M=2,Q' > | 2 
£( J, M=l, 0=1) - E(J>, M=2, fi 7 ) - frwj. 




E 


/ ^ 12 


| < J,M=-l,H=-l|d + |J',M=-2,H' > 

E( J, M= — 1, Sl= — 1) — E( J', M= - 2, IT) + ~~ 

I < j,m=- i , n =- i | d _| j',m=o,w > | 2 

“ E(J, M= — 1, ri= — 1) — E(J', M= 0,12') - ?iw_L ' 

( 3 ) 


SAE 1 =£l/4x 

| < J,M=l,Q=-l\d + \J',M = 0,fi' > | 2 
f-' E{ J,M=l, 12= - 1) - £( J', M= 0,IT) + 


(E 


j',n 


E 


< j, m=i, ci= — i|d_|j', m= 2,ry > I 2 

f-' E{ J, M=l, 12= - 1) - A( J', M= 2, IT) - hujj_ 


of magnitude whereas Sd e (sys ) is almost unaffected by 
this interaction. See also influence of perturbation by 
J — 2 level on J = 1 g-factors in Refs. SEMI [El- 

Tables m and nu list the calculated Sd e (sys) and 
Sd e (sys) as a functions of u>± and £ for J = 1 and J = 2, 
correspondingly. Though for smaller £ the £± value will 
be smaller as well, for the calculation I take the same 
£± = lOmV/cm given in Ref. [II for all £. Using the 
fact that Sd e (sys) and Sd e (sys) are quadratic functions 
of £± the results can be easily recalculated for any £±. 
For static magnetic field the value B = 40mG used in the 
experiment [5J is used. One can see that Sd e (sys) is two 
orders of magnitude larger for J =2 than for J=1 though 
much smaller than the current limit on d e . 

Calculation for w_l/27t less than 250 kHz is not per¬ 
formed due to the limited computational accuracy. For 
smaller lo± one can expect further decreasing of Sd e (sys) 
and 5d e (sys). Each term in Eqs. (13141) has form 

Retaining terms up to the third 


»(0 


»(0 


o u (I)+fi«l a u (i)-hu>±‘ 

order in w_l we have 

SAE u( - l) k, -2- S “ (i) 


hui± 


j',n 


A. 


u{l) 


A. 


u(l) 


n 3 co 3 j_ 

A 40 


( 5 ) 


-E 

j\n 


\ < J,M= 
, E(J, M= - 


-E 

jqn' 


| < J, M 
E(J, M= - 


---l,Sl=l\d+\J , ,M=-2,SV > | 2 
1, H=l) - E(J', M= - 2, IT) + huj ± 

= - l,Q=l\d-\J',M=0,Cl' > | 2 
-1, £1=1) - E(J\ M= 0, £1') - hu ± ’’ 


( 4 ) 


Major contribution to SAE U ^ comes from coupling of 
states with the same J. The most simple is the pic¬ 
ture for J= 1 state. |<7=1, M=l, £1, n > interact with 
| J=l, M= 0, £1, n +1 > and | J = 1, M= — 1, f l,n > with 
| J, M= 0, £1, n— 1 >. Note, that H mo i can only couple the 
states with the same n, whereas Hint couples the states 
with An = ±1. Energies of states | J, M=0, £1, n+1 > 
and | J, M=0, £1, n— 1 > are different by and 

states \J=1, M=l,£l,n > and | J=l, M= — 1, — fi, n > 
by 2 [iB. This leads to different energy denominators 
in eqs. m and results in different energy shift for 
\J=1, M=l,£l,n > and |J=1,M= — 1, —f l,n >. How¬ 
ever, for J= 1 level, it was shown in Ref. [H} that (consid¬ 
ering the interaction with \J=l, M=0,£l > states only) 
provided tensor Stark (A Est = E(J,M= ±1,H) - E(J, 
M=0,f2)) and Zeeman splitting are the same for upper 
and lower component of the H-doublet the AC Stark 
shifts 6AE U and SAE 1 will also be equal. This allows 
one to reject systematic errors due to geometric phases 
by performing measurements in both £1— doublet states. 
However the tensor Stark splittings do not coincide ex¬ 
actly. Also including interaction with other states lift the 
degeneracy. It is particularly important to include the in¬ 
teraction with the neighbor rotational levels. The latter 
interaction increases value for Sd e {sys ) on several orders 


Formulae for B u m and A u ^ for J= 1 are given below. Eq. 
© explains the fact that SAE U T decreases linerly with 
small uj_ i listed in Tables [Q and HU Similarly to 5AE U ^ 
the major contribution to difference SAE 1 — 5AE U comes 
from coupling of states with the same J (terms with 
J' = J in Eqs. EBP )- However, important role plays 
the perturbation by the closest rotational levels which 
makes matrix elements and denominators for upper and 
lower components of Sl-doublet slightly different. Let us 
consider this effect for the simplest case the J— 1 level. 
Without perturbation by the J = 2 level the parameters 

A = AE S t =-£ < J=l,M=l,£l\d z \J=l,M = l,£l> 

= -£dM£l/J(J+ 1) 

( 6 ) 

and 

B = -£j_/2 < J=l,M=l,£l\d + \J=l,M'=0,£l> 

£ ± dn v / (J-M+l)(J+M) 

2 J(J+ 1) ( ’ 

up to the sign are the same for upper and lower H-doublet 
levels. A Est is positive for upper and negative for lower 
levels. Note that dipole moment d < 0. Eqs. EEP ex- 
plain the fact that SAE U ^ decreases quadratically with 
£. 

Perturbation by J = 2 changes the parameters: 

A U{1) = A E st + S'A&gp + S 2 AE^\ (8) 


B u (i) = B + S 1 B u (j) + 5 2 B u (i). 


(9) 
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TABLE I. The Sd e (sys) (in units 10 29 e-cm ) and 5d e {sys) TABLE II. The 5d e (sys) (in units 10 29 e-cm ) and Sd e (sys) 
(in units 10 -34 e-cm ) calculated for the J = 1 H 3 Ai state in (in units 10 _34 e-cm ) calculated for the J = 2 H 3 Ai state in 

232r-pj^l6Q 232^^160 


£ = 4MHz wx/ 27r = 

(V/cm) 8d e (sys) 8d e (sys) 8d e (sys) 

= 1MHz cjx/27t = 250kHz 
Sd e (sys) 8d e (sys) 8d e (sys) 

£ 

(V/cm) 

Wx/27r = 4MHz 
8d e (sys) 8d e (sys) 

Wx/27r = 1MHz 
8d e (sys) 5d e (sys) 

CJx/27T = 250kHz 
5d e (sys) 8d e (sys) 

20. 

-1299. 

-2144. 

-313. 

-31. 

-78. 

-0.48 

20. 

-1996. 

-213498. 

-320. 

-1663. 

-78. 

-33. 

30. 

-565. 

-608. 

-139. 

-9.1 

-35. 

-0.14 

30. 

-658. 

-42436. 

-140. 

-495. 

-35. 

-9.9 

40. 

-315. 

-253. 

-78. 

-3.8 

-19. 

-0.059 

40. 

-342. 

-15657. 

-79. 

-211. 

-20. 

-4.2 

50. 

-201. 

-128. 

-50. 

-1.9 

-12. 

-0.030 

50. 

-212. 

-7546. 

-50. 

-107. 

-13. 

-2.1 

60. 

-139. 

-74. 

-35. 

-1.1 

-8.6 

-0.018 

60. 

-144. 

-4229. 

-35. 

-61. 

-8.7 

-1.2 

70. 

-102. 

-47. 

-25. 

-0.70 

-6.3 

-0.011 

70. 

-105. 

-2613. 

-26. 

-39. 

-6.4 

-0.78 

80. 

-78. 

-31. 

-19. 

-0.47 

-4.8 

-0.0073 

80. 

-80. 

-1728. 

-20. 

-26. 

-4.9 

-0.52 

90. 

-62. 

-22. 

-15. 

-0.34 

-3.8 

-0.0053 

90. 

-63. 

-1204. 

-15. 

-18. 

-3.9 

-0.37 

100. 

-50. 

-16. 

-12. 

-0.24 

-3.1 

-0.0038 

100. 

-51. 

-872. 

-13. 

-14. 

-3.1 

-0.27 

110. 

-41. 

-12. 

-10. 

-0.18 

-2.6 

-0.0028 

110. 

-42. 

-653. 

-10. 

-10. 

-2.6 

-0.20 

120. 

-35. 

-9.2 

-8.6 

-0.14 

-2.2 

-0.0022 

120. 

-35. 

-501. 

-8.7 

-7.7 

-2.2 

-0.15 

130. 

-30. 

-7.3 

-7.3 

-0.11 

-1.8 

-0.0017 

130. 

-30. 

-393. 

-7.4 

-6.1 

-1.9 

-0.12 

140. 

-26. 

-5.8 

-6.3 

-0.088 

-1.6 

-0.0014 

140. 

-26. 

-314. 

-6.4 

-4.9 

-1.6 

-0.098 

150. 

-22. 

-4.7 

-5.5 

-0.071 

-1.4 

-0.0011 

150. 

-22. 

-255. 

-5.6 

-4.0 

-1.4 

-0.080 

160. 

-20. 

-3.9 

-4.9 

-0.059 

-1.2 

-0.00093 

160. 

-20. 

-210. 

-4.9 

-3.1 

-1.2 

-0.063 

170. 

-17. 

-3.2 

-4.3 

-0.050 

-1.1 

-0.00078 

170. 

-17. 

-175. 

-4.3 

-2.7 

-1.1 

-0.055 

180. 

-15. 

-2.7 

-3.8 

-0.043 

-0.96 

-0.00067 

180. 

-15. 

-147. 

-3.9 

-2.2 

-0.97 

-0.044 

190. 

-14. 

-2.3 

-3.4 

-0.034 

-0.86 

-0.00053 

190. 

-14. 

-125. 

-3.5 

-1.9 

-0.87 

-0.038 

200. 

-12. 

-2.0 

-3.1 

-0.030 

-0.77 

-0.00047 

200. 

-13. 

-107. 

-3.1 

-1.6 

-0.78 

-0.032 


S^AEst is the correction to A Est due to shift¬ 
ing down |J=1,M=± 1,11 > (|J=1,M=0,12 >) lev¬ 
els when interacting with J= 2. 5 1 AE'^) is negative. 

It decrease (increase absolute value) A Est for upper 
(lower) 12-doublet levels. In turn 6 2 AEg%) is posi¬ 
tive. It increases (decreases absolute value) A Est for 
upper (lower) 17-doublet levels. S 1 ^B is the correc¬ 
tion to B due to the perturbation of the wavefunction 
| J=l, M= ± 1,12 > (| J=l, M= 0, 12 >) by | J=2, M= ± 
1,17 > (|J=2, M=0,12 >) one. It is shown in AP¬ 
PENDIX that for J=1 level the corrections 5 1 AEg^), 
S 2 AEgip , 5 1 B U p), and S 2 B U p) are correlated in such a 
way that 


r>2 

U u(l) 

A 2 

u(l) 


B 2 

A 2 ' 


( 10 ) 


Eq. (HOD is correct up to the second order in small param¬ 
eter A Egr/AE rot , where A E rot = E(J= 1) — E(J= 2) is 
energy difference between the first and second rotational 
levels. Due to Eq. m the linear term in the difference 
Sde(sys) = SAE 4 £^ EU is canceled and Sd e (sys), in the 
leading order, is a cubic function of for J= 1. This 
behavior can be seen from the data in Table |H Depen¬ 
dence of the Sd e (sys) for J =2 level on uj± has also nearly 
the cubic character. 

The calculations confirm that the experiment on ThO 
H 3 Ai state is very robust against systematic errors re¬ 
lated to geometric phases. Developed code can be applied 
for calculation of molecules in an ion trap at presence of 
rotating field QQ. 


The author is grateful to A. V. Titov for useful dis¬ 
cussions. Codes for calculations of diatomic molecules 
were developed with the support of the Russian Science 
Foundation grant (project No. 14-31-00022). Calcula¬ 
tion of the ThO molecule were performed with the sup¬ 
port of Saint Petersburg State University, research grant 
0.38.652.2013 and RFBR Grant No. 13-02-0140. 


APPENDIX 


In the first order in the small parameter 
AE ST /AE rot ~ £d/AE rot for S 1 AE^\ S 2 AE^\ 
5 1 B U V\ and S 2 B U ^ we have 

<5 4 A E$p = 

£ 2 /AE Iot | < J=l,M=l,n\d z \J'=2,M = 1,17 > | 2 = 
£ 2 d 2 ((J+l) 2 - M 2 )((J+1) 2 - 12 2 ) 

A E rot (2J+1)(2J+3)(J+1) 2 

S 2 AEg { rl) = 

-£ 2 /AE mt \ < J=l,M=Q,n\d x \J'=2,M = 0,17 > | 2 = 

£ 2 d 2 (J+l) 2 ((J+l) 2 — 12 2 ) 
_ AE rot (2J+1)(2J+3)(J+1) 2 ’ 


S 1 B u(i) = £±/2 < J'=2,M=l,£l\d+\J=l,M'=0,Q > x 
£/AE mt < J=1,M = l,Ct\d z \J'=2,M = 1,17 >= 

£±£d 2 /(J+M)(J+M+1)((J+1) 2 -17 2 ) w 
_ 2AE rot y (2J+1)(2J+3)(J+1) 2 X 

/ ((J+l) 2 — M 2 )(( J+l) 2 — 17 2 ) 
V (2 J+l)(2J+3)( J+l) 2 
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S 2 B u(l) =£_!_/2 < J=l,M=l,n\d+\J'=2,M'=0,n> x 
£/AE Iot < J=l,M'=0,n\d z \J'=2,M'=0 >= 

£±£d 2 /(J-M+2)(J-M+1)((J+1) 2 - fl 2 ) 


2AE r , 


(2J+1)(2J+3)(J+1) 2 


(J+ 1 ) 2 ((J+ 1 ) 2 - n 2 ) 


(2J+1)(2J+3)(J+1) 


2 ' 


therefore 6^AE{^p<0 and 


AE rot is negative, 

S^AEg^>0. 

Then retaining terms up to the first order in 
AEst / AE rot we have 


B 2 ,,, 

A 2 

«(0 


(i? + S 1 B u (i' ) + 5 2 B u (i )) 2 
{AEst + 5 1 AE‘^} + 8 2 AE'$) 2 
B 2 + 2B6 l B u(l) + 2BS 2 B u{1) 


AEg T 


2AE st 6 1 AE 


*(0 

ST 


2AE st 5 2 AE 


*(0 

ST 


77(7+1)?- 


£ 2 d 2 M 2 i 2 2 

(J(J+1)) 2 


(i + [(j+M+i) - (J+i)</^gp] A-(J)/n) 


B 2 (l+ 


(i +Afif(j)/n) 

(J+M+l) — (J+l)y/ 


J-M +2 

J+M 




A 2 (l + M/\(J)/ff) 


( 11 ) 


where 


K (J) = 2 


£d J(J+1)((J+1) 2 -ft 2 ) 
AE Iot (2J+1)(2J+3)(J+1) 2 ' 


Substituting M= 1, J=1 to Eq. dill we have got Eq. 
(flOl) . Eq. (fill) is obtained for M=+l level. The re¬ 
sult is the same for M=— 1. Taking into account that 
| J=l, M= + 1, n= + 1 > (I J=l, M= + l i n=-l>) cor¬ 
responds to the upper (lower) fl-doublet level, for the 
next term we have 

B 2 U _ B 2 (l + K) _ B 2 { 1 + 5 ^) 

A4 ^4 {1 + 2K) A4 (1+2 ^-)’ 

Bf B 2 (l — K) B 2 {l-^- t ) 

Af A 4 (l — 2K) A\l-2^± 7t y 

As an example, below are the parameters calculated 
for £ = 110 V/cm. 


2 hBl 
4 £ eS A 2 u 


= 10.33057851 x 10~ 29 


e-cm 
2ttMHz ’ 


2 TiB 2 

4^4f 


= 10.33057851 x 10" 29 


e-cm 
2ttMHz ’ 


2 h 3 B 2 u 
4£ ef fAi 


0.00080156 x 10“ 29 


e-cm 

(2ttMHz) 3 ’ 


2h 3 B 2 
4 £esAf 


0.00080343 x 10" 29 


e-cm 

(2ttMHz) 3 ' 
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